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Detection and correction of phase jumps in a phase sequence _ 



Field of the invention 

The invention relates to a communication system comprising at least a 
transmitter and a receiver intended to receive symbols coming from a PSK modulation, and 
comprising estimation means for estimating a frequency error relating to a symbol based on a 
sequence of symbol phases. The invention also relates to a receiver intended to be used in 
such a communication system. 

The invention also relates to a method of estimating a frequency error relating 
to a received symbol, coming from a PSK modulation, based on a sequence of symbol 
phases, and to a method of detecting and correcting phase jumps in an initial phase sequence 
of symbols coming from a PSK modulation. 

The invention finally relates to computer programs comprising instructions for 
implementing such methods. 

Background of the invention 

Such an estimation algorithm of a frequency error relating to a received 
symbol based on a phase sequence is described, for example, in paragraph 4.1 (page 107) of 
the article «Feedforward Frequency Estimation for PSK : a Tutorial Review» by 
M. Morelli and U. Mengali, published in the journal «European Transactions on 
Telecommunications, vol. 9, no. 2, March - April 1998». This algorithm is known by the 
name of Tretter algorithm, or least squares method. 

To obtain such a phase sequence, it is known that a phase estimation algorithm 
is used, which estimates the phase relating to a received symbol on the basis of decisions 
made on various received symbols. For example, the expectation maximization algorithm is 
used, which is described in paragraph 3.3 of the conference report of the -international 
Conference on Communications, New Orleans, USA, 1-5 May, 1994, vol. 2, pp. 940 and 
945», entitled «Comparison between digital recovery techniques in the presence of 
frequency shift» by F. Daffara and J. Lamour. 
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The problem posed is the following: in a PSK modulation comprising 2 n 
points, two adjacent points have a phase difference of . When the frequency error is such 
that, based on a certain symbol, an error is made in the decision, this error is translated by a 
phase jump of +-^- in the sequence of phase estimates obtained. The number of phase 

jumps that may be obtained depends on the number of symbols contained in the packet and 
the initial frequency difference. 

When the sequence of phases, which are used for applying the Tretter 
algorithm, includes one or various phase jumps, the frequency estimate obtained is 
inaccurate. The invention notably has for its object to provide a solution to this problem. 

Summary of the invention 

For this purpose, a communication system according to the invention and as 
described in the opening paragraph is characterized in that said receiver comprises 
calculation means for calculating a phase sequence, called initial sequence, based on 
decisions made on symbols, and means for detecting and correcting phase jumps in this 
initial sequence, to supply a phase sequence, called final sequence, to said frequency error 
estimation means. 

Advantageously, said means for detecting and correcting phase jumps 

comprise: 

- modifying means for modifying said initial sequence, so as to produce a plurality of 
modified sequences which each compensate for a phase jump configuration, 

- calculation means for calculating straight line equations which determine the initial 
sequence and the modified sequences, 

- calculation means for calculating for the initial sequence and the modified sequences a 
mean difference between the initial or modified phases and the phases produced by the 
corresponding straight line equation, said final sequence being formed by the sequence 
whose mean difference is minimal. 

Brief description of the drawings 

These and other aspects of the invention are apparent from and will be 
elucidated, by way of non-limitative example, with reference to the embodiment(s) described 
hereinafter. 
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In the drawings: 

- Fig. 1 is a diagram of an example of a communication system according to 

the invention, 

5 - Fig. 2 is a flow chart describing the steps of a phase estimation method, 

- Fig. 3 is a flow chart diagram describing the operations used by the phase 
jump correction and detection means according to the invention, 

- Fig. 4 is a representation of a phase sequence including a phase jump, 

- Fig. 5 is a representation of a phase sequence including two phase jumps, 
1 0 - Fig. 6 is a representation in the form of curves of the results obtained with 

the invention. 

M> 

q Description of preferred embodiments 

.-- 

In Fig. 1 is represented an example of a communication system according to 
%0 15 the invention. This communication system comprises interactive user terminals 1, which are 
jr transmitters within the meaning of the invention, and a head-end station 2, which is a receiver 

within the meaning of the invention. The head-end station 2 transmits signals in a first 
O frequency band Ku (12-14 GHz). These signals are relayed to interactive user terminals 1 by 

ry 

m a satellite 3 . The interactive user terminals transmit signals in a second frequency band Ka 

20 (20-30 GHz). These signals are relayed to the head-end station 2 by the satellite 3. 

f|f Each terminal 1 comprises a data source 10 and channel coding means 12. The 

channel coding means deliver packets of N symbols which contain preamble symbols and 
data symbols. These packets are then transmitted to filter means 13 and, finally, to 
modulation means 14, which use a local oscillator that has a frequency f c . 
25 The head-end station 2 comprises demodulation means 20 which use a local 

oscillator that has a frequency f c +Afo/Ts (where Afo is a normalized frequency difference 
relative to the symbol frequency, and Ts is the duration of the symbols), and an initial phase 
@o- The head-end station 2 also comprises filter means 21 and sampling means 22, which 
sample the output signal of the filter 21, to deliver symbols called received symbols. The 
30 received symbols are transmitted to frequency recovery means 23 which estimate the 

normalized frequency difference Afo and correct the received symbols to compensate for the 
estimated difference A/o • The frequency recovery means 23 deliver frequency-corrected 
symbols. These frequency-corrected symbols are transmitted to phase recovery means 24. 
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The phase-correeted and frequency-corrected symbols are finally transmitted to channel 
decoding means 25 which deliver data 27. 

In the example described here, the phase recovery means 24 comprise phase 
and frequency estimation means 50 and phase correction means 5 1 . The phase estimation 
means 50 are described with reference to Fig. 2. They are formed by a loop intended to be 
passed through L times. In the following of the description the index m (m=l to L) is a loop 
counter. Each value of m thus corresponds to one loop path. The loop comprises: 

- conventional phase estimation means 52 (for example, expectation maximization) for 
producing an initial phase sequence SI relating to a symbol packet r£ (TO_1) (m=l to L, and 
k=l to q, where q<N), 

- means 53 for calculating a frequency error relating to said packet, based on the initial 
sequence SI, 

- means 54 for correcting the frequency of the symbols of the packet, to correct said 
frequency error, 

- loop means 55 which supply the frequency-corrected symbols r£ (m) to phase estimation 
means 52 for a following path through the loop, 

- and loop output means 56 which supply to the phase correction means 5 1 the phases 
estimated by the phase estimation means 52 (or, directly, the correction to be made), and 
the symbols to be corrected. 

The calculation means 53 calculate a frequency error estimate Af m relating to 
the symbols r/t (m_1) of a same packet, based on the initial phase sequence SI produced by the 
phase estimation means 52. The values of this phase sequence SI are between -co and +oo. 
The calculation means 53 comprise: 

- means 60 for detecting and correcting phase jumps to correct this initial sequence SI and to 
deliver a final sequence S2, 

- means 62 for calculating the slope of a straight line which is closest possible to the values of 
the sequence S2, by applying the Tretter algorithm. The slope obtained forms a frequency 
error estimate Af m relating to the symbols of the packet. It is this frequency error that is 
transmitted to the frequency correction means 54. The symbols obtained after frequency 
correction, rk (m) = r^" 1 '^ .e~ 2 *^ k ~ hfm , are transmitted to the phase estimation means 52 for a 
new path through the loop. During the last path through the loop (m=L), the symbols r^ L ^ 
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to be corrected and the phase cprrection e J ® k to be made to these symbols are 
transmitted to the phase correction means 51 . 

The means 60 for detecting and correcting phase jumps are represented in Fig. 
3. They comprise: 

5 - means 1 00 for modifying the initial sequence to compensate for a plurality of phase jump 
configurations; the correction means 100 produce a plurality of modified sequences which 
correspond each to the correction of a phase jump configuration; 

- means 1 10 for calculating straight line equations which determine the initial sequence and 
the modified sequences, 

10 - calculation means 120 for calculating for the initial sequence and the modified sequences a 
mean difference between the initial or modified phases and the phases produced by the 
corresponding straight line equation, said final sequence being formed by the sequence 
I whose mean difference is minimal. 

n 

pi The use of the means 60 for detecting and correcting phase jumps is different 

•ft 1 5 and depends on the number of phase jumps one wishes to correct. But the method applied 

=P remains the same. The invention can thus be applied to any number of phase jumps. Now two 

M 

= examples of embodiment will be described of the means 60 for detecting and correcting 

h 

-1] phase jumps: 

- a first example, in which the means 60 for detecting and correcting phase jumps are 

a 

. 20 intended to correct a single phase jump of ±2- ; 
ft! 2 

- a second example, in which the means 60 for detecting and correcting phase jumps are 

intended to correct two phase jumps of ±-| in the same direction; this second example 

corresponds to the most probable case where the transmitted packets are ATM cells of 53 
octets. 

25 In the examples that will be described, the modification means 1 00 modify the 

initial sequence SI phase by phase. But, for diminishing the number of calculations to be 
performed, it is possible to modify the sequence phase-group by phase-group. This is 
equivalent to compensating only for certain configurations of phase jumps. 

The first example of embodiment of the means 60 for detecting and correcting 

30 phase jumps is described with reference to Fig. 4. In Fig. 4 is represented an example of an 
initial sequence S 1 comprising a phase jump of +-|- . The means 60 for detecting and 
correcting phase jumps have for their function to detect the position and the direction of the 
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phase jump and then correct it. Therefore, as indicated in Fig. 3, they perform the following 
operations: 

a) The Tretter algorithm is applied to the initial sequence SI formed by phases cpj (j=0 to q-1) 
to obtain the straight line equation D 0 which determines this sequence. This equation is 
5 written as: 

D 0 : y=ao.x+b 0 with: 

a Q =a.S'-j3S and b Q =y.S-J3S' 

where S = §^ , S'=^j^j ,« = 12/(q.(q 2 -1)) , P = 6/(q.(q + l)) and 

j=o j=o 

r = 2(2q-l)/(q(q + l)) 

1 0 These expressions can easily be derived from the calculations shown on pages 523 and 524 
M of the title «Numerical Recipes in C, the art of scientific computing, second edition» by 

p W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery, published by Cambridge 
5 University Press in 1 995, while considering that the uncertainty of the phases is constant 
SI whatever j. 

P 15 b) An initial mean difference (cr 0 ) 2 is calculated between the phases q> j of the initial 
sequence SI and the phases y(j) coming from the straight line equation D 0 . . 




c) The initial sequence SI is run through point by point by starting from the end (symbol of 
rank q-1); the index i is a counter that indicates the position of the phase jump (i=q-l, ...,0). 

20 d) With each step the phases q>\ (j = q - i, . . .,q - 1) are modified by + ^ , so that a modified 
sequence Q + is obtained. This sequence Q + is thus constituted by phases 
f <p. pour j = 0,...,q-i-l 

^ j <p \ +f-P° ur J = vq-i 

e) With each step one straight line equation Dj + is calculated, which straight lines determine 
the modified sequence Q + . This equation is written as: 
25 Di 4 : y=(a,) + .x + (b;) + with : 
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aj + =«.(S j ') + -/?.(Sj )+ and b s + =r.(Sj) + -A(Si') + 

where (S,) + = = S + I.± and (S,') + = ][>/ -S4-§J(q-^) 

j=0 z j=o 

that is to say, 

a;=a o +a.|i.(q-^)-Af = a 0 +(A;) with(A + ) = a.|i.(q-^)-^.| 
and b i + =b o +yi.|-A|i.(q-i|i) = b o +(B i + ) with (B i + ) = y i.|-)?.|i(q-^) 
f) For each modified sequence Q + is calculated a mean difference (cr; + ) 2 between the phases 
97 j + of the modified sequence Q + and the phases y(j) coming from the straight line equation 



{m + ) 2 =-Y\<Ps + -[(ai) + -i + (bi) + ] 
g) The operations c) to f) are repeated while the phases (j = q - i,...,q - 1) of the initial 
sequence of - are modified. For each value of i, another modified sequence Cf is 



obtained. It is formed by the phases (p = 



for j = 0,...,q-i-l 

j 

< p.-^forj = g-i,...q-l 



h) The final sequence S2 is formed by the sequence whose mean difference is minimal. 

For a less complex use, the mean differences (cr*) 2 are calculated based on 
the initial mean difference (a Q ) 2 . One obtains: 

(a+) 2 =a Q 2 +2x <p> -2.(Bj + ).S-2(Ai + ).S' 

j = q-i 

+ b 0 [(Aj + ).q.(q-l) + 2.(Bi + ).q-^i] 

+ a 0 [ (Ai+) - q - (q 3 1) - (2q " 1) + (B i + )-q-Cq-l)-f.i.(2q-i-l)] 

The mean difference (cr j ~) 2 is obtained by replacing in the expression of (er* ) 2 : 
xby-n, (A") by - (Af) and (B ; + ) by - (B ; ") . 
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The second example of embodiment of the means 60 for detecting and 
correcting phase jumps is described with reference to Fig. 5. In Fig. 5 is shown an example of 
an initial sequence SI comprising two phase jumps of + y .At step d) the phases are 

modified by +|- for the p symbols from rank q-k-p to q-k-1, and by + n for the k symbols 

from rank q-k to q (k varies between 1 and q and p varies between 1 and q-k). 
The modified sequences obtained at step d) are thus written as: 
<p . for j = 0,...,q-k-p-l 
<p } +yfor j = q-k-p,...,q-k-l 
<p. +7t for j = q-k,...,q 

(p. for j = 0,...,q-k-p-l 
cp -y for j = q-k-p,..., q-k-1 
q> -7i for j = q-k,...,q 
And the straight line equations calculated in step e) are written as: 
D p ,k + : y = (a P ,k + ).x + (b p ,k + ) with: 

( a P/k ) + =*-( S p,k'> + -^ S p,k ) + and ( b p,k> + = ^ S p,k )+ -?< S P P + 
where (S pk ) + =X^j + = S + p.y + lor and 

(S pi 'r=|j^;-S' + f.p.(q-k-^) + ^.(q-^) 
that is to say, 

a p ;=a o+ a[|. P .(q-k-^) + ^.k.( q -^)]-^[p| + k^] = a o+ (A pk + ) 

with (A + ) = a£ .p.(q -k-2±±) + *Jc(q - ^±1)] - fi\p % + k*] 
p,k 2 2 2 2 

andb * = b +r[p^ + k^]-^-P-(q-k- £ ^-) + ^.k.(q-^)] = b +(B + ) 

P ,k o 2 2 2 2 u p,K 

with(B + ) = r[p^ + k^-^.p.(q-k-^) + ^.k.(q-^)] 
p,k 2 2 2 Z 

For a less complex implementation, the mean differences 0 P/k + ) 2 are calculated 
based on the initial mean difference (cr ) 2 . The following expression is obtained: 
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Kk + ) 2 =cr 0 2 ^ S ?j '"l^j -2.(B p p.S-2(Ap, k + ).S' 
j=q-k j=q-k-p 
+ b 0 [(A p ,k + ).q.(q - 1) + 2.(B pk + ).q - «(p + 2k)] 

(A P> k + ).q.(q-l).(2q-l) +(Bpk + )q(q _ 1) _| p(2q _ 2k _ 

+ q.B p/ ) 2+ (A P , k+ ).(B p/ ).q.(q-l) ^ A ^> 2 ^ ( r 1) - (2q - 1) 
- 1 -(Ap,k + ).p.(2q - 2k - p - 1) - ^.(A p ,k + ).k.(2q - k - 1) - *.(B pk + ).(p + 2k) + ^-.(p + 4k) 
The mean difference (cr p k ~) 2 is derived from the expression (o" P/k + ) 2 b Y replacing 

n by -w, (A p /) by - (A pk ") and (B pk + ) by - (B pk ") . 

It is evident that the means that have just been described are calculation means 
Q 5 advantageously used in the form of a computer program intended to be executed by a 
H microprocessor placed in the receiver. 

^0 In Fig. 6 are shown the results obtained thanks to the invention in a system as 

J described with reference to Figs. 1 and 2. The curves of Fig. 6 represent the packet error rate 

(PER) plotted against the signal-to-noise ratio (SNR) in the following three cases: 
• '10 - curve R3 : when the phase jumps are not detected; 

ry - curve R2: when the phase jumps are detected and corrected with the method according to 
the invention; 

ry - curve Rl : for a Gaussian channel (no imperfection as regards either phase or frequency). 

These curves have been obtained via: 
15 - a QPSK modulation, 

- a random initial frequency error Af 0 comprised between -1% and +1% of the symbol 
frequency, 

- a random initial phase error 9 0 comprised between - % et + n , 

- packets which comprise a known preamble of 48 symbols [ai, a4g] and a payload part 
20 [&49, a 592 ] of 544 symbols, 

- a frequency recovery algorithm (utilized by the means 23), which uses the packet 
preambles, and which is applied before the phase estimation, so that the residual frequency 
error Afi before the phase estimation is lower than or equal to 0.3% of the symbol frequency, 

- phase estimation means 50 which are formed by a loop which is passed through L=2 times, 
25 - on transmission, a Reed Solomon coding and a convolution coding which is punctured in 64 

states. 
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The invention is not restricted to the embodiments that have just been 
described by way of example. More particularly, it relates to any type of receiver that utilizes 
a frequency error correction device of the Tretter type and, upstream of this device, a phase 
estimation device based on decisions made on the received symbols. 

Moreover, the number of phase jumps that may be detected and corrected by 
applying this method is arbitrary. When the transmitted data packets are longer, it is possible 
to have more than two phase jumps. In that case, in order not to complicate the calculations 
too much, one advantageously chooses to divide the data packets into smaller portiona so as 
to be in the same position again where one has a maximum of two phase jumps per packet 
portion. The method described above is thus applied to each portion of the packet. It provides 
the position and the direction of the various phase jumps. The phases are then corrected by a 
multiple of y as a function of the direction of the jumps and their number. Then, the Tretter 

algorithm is applied to this corrected sequence to obtain the final frequency estimate. In order 
to avoid a jump occurring between two portions of a packet, it is desirable to provide an 
overlap between the various portions of the same packet. 

The method proposed is generally used for any MPSK modulation by 
considering phase jumps of ±2tt/M (instead of + a 1 2 for a QPSK modulation). 



